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ABSTRACT 

A dual beam apparatus was developed which simultaneously 
measured particle size (0^2) the entrance and exit of an 
exhaust nozzle of a small solid propellant rocket motor. 

The diameters were determined using measurements of diffrac- 
tively scattered laser power spectra. The apparatus was 
calibrated by using spherical glass beads and aluminum oxide 
powder. Measurements were successfully made at both loca- 
tions. Because of the presence of char agglomerates in the 
exhaust, continued effort is required to improve the grain 
design in order to obtain consistent "across nozzle" data. 
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I. 



INTRODUCTION 



Aluminum is used as a fuel additive in solid propellants 
to increase the specific impulse of the rocket and to sup- 
press high frequency combustion instability. Although there 
is an increase in specific impulse, the specific impulse ef- 
ficiency may actually decrease when compared with the effi- 
ciency of the same base propellant without aluminum [Ref. 1] . 
This performance loss has been largely attributed to incom- 
plete combustion of the aluminum and/or the formation of 
aluminum oxide particles in the motor cavity which leads to 
two-phase flow losses (i.e. velocity and thermal lags between 
the particles and the gas) . The magnitude of the losses as- 
sociated with these mechanisms has been reported to be 1% I 

sp 

loss for 10% unburned aluminum and 2-10% loss in I for two- 

sp 

phase flow effects [Ref. 2]. These losses are predicted 
analytically and at present there exists no adequate theo- 
retical model relating particle size to propellant and motor 
parameters. Therefore, the accuracy of performance predic- 
tion is very dependent upon the accuracy of particle size 
data obtained experimentally. 

The process by which aluminized fuel additives proceed 
to final combustion products is exceedingly complex. The 
final particle size is a function of several factors: ori- 

ginal particle size, propellant properties, operating 
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environment (pressure, etc.) and the nozzle design and throat 
size [Ref . 2 ] . 

Kincaid and Derr [Ref. 3] have presented a comprehensive 
description of the process. A summary of their explanation 
is presented here. They divide the rocket motor into five 
zones: (1) at or near the propellant surface (0.5 to 5 mm) ; 

(2) within the motor cavity; (3) within the converging sec- 
tion of the nozzle; (4) at the nozzle throat; (5) within 
the diverging section of the nozzle. In zone (1) , each par- 
ticle is heated as it emerges at the surface. Its metal core 
melts, expands and eventually ruptures its surrounding oxide 
layer. Some of the original metal particles may ignite and 
leave the burning surface immediately. Other particles may 
accumulate on the surface in groups of 100 to greater than 
1000 — sometimes a million [Ref. 4]. These groups leave the 
surface either as clumps or burning droplets, called agglom- 
erates. The aluminum oxide particles produced by aluminum 
combustion process have a bimodal size distribution. Fine 
material (1 ym) results from vapor phase combustion. This 
size is generally considered to be invariant and upon com- 
pletion of aluminum combustion amounts to approximately 80% 
by weight of the total condensed phase products. Coarse 
material (100 ym) arises from breakup of the residual cap of 
AI 2 O 2 on each burning particle. This large cap remains after 
the combustion of the metal has completed. The size of these 
particles is dependent on; propellant composition, pressure. 
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flow environment and the original size of the particle or 
agglomerate. It is believed that the degree of agglomera- 
tion on the surface will increase for highly aluminized pro- 
pellants. Hence the proportion of coarse particles would 
increase . 

Some oxidation occurs in zone (2) , but the precise mech- 
anism of oxidation is not known. In addition, some growth 
of AI 2 O 2 particles by condensation and coagulation occurs in 
this zone. Behavior within zone (2) has been inferred from 
laboratory data. It has been suggested that the behavior in 
zone (2) would be dependent mainly on the size of the parti- 
cles entering the zone. With increased agglomeration in zone 
(1) and short residence time in zone (2) the particles could 
enter the nozzle and be ejected prior to complete combustion. 
Another complication could be settling of the larger particles. 

In zones (3) , (4) and (5) , size is affected by wall col- 

lisions, growth of particles, breakup of particles, condensa- 
tion and solidification. Zone (3) is a region of accelerating 
gas velocity. When the two-phase flow of zone (2) enters 
this zone, collisions occur which lead to the growth of 
larger particles. As these particles accelerate, their 
break-up is also expected. Because of the change in cross- 
sectional area, the flow must turn and wall collisions occur. 

As the particles enter the throat of the nozzle (zone 4) , 
particle breakup has been observed [Ref. 5] . As the flow 
expands through zone (5) , it is cooled, and condensation and 
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solidification of the alumina occurs. The problem of two- 
phase flew losses occurs within zones (3) -(5) when the con- 
densed phase material cannot maintain both velocity and 
thermal equilibrium with the gaseous phase. According to 
theory, condensed material smaller than 1 ym in diameter en- 
tering the nozzle will produce minimal two phase flow losses 
even for highly aluminized solid propellants. 

At this time the Air Force Improved Solid Performance 
Program (SPP) uses a particle size model developed by Hermsen 
[Ref. 6], based on ^ weight (or DeBroucker) mean diame- 

ter. This model is based on correlation of a large amount of 
empirical data from many sources, most of which employed par- 
ticle analysis of collected exhausts. Particle break-up and 
fragmentation mechanisms within the nozzle utilize a model 
based on the critical Weber number concept [Ref. 6] . How- 
ever, test cases have indicated that assumptions made in the 
model yield unrealistic particle distributions [Ref. 7]. 

Very little experimental information is available that re- 
lates nozzle geometry and inlet particle size to the behavior 
of the particles as they pass through the nozzle. Light 
scattering techniques were not considered in obtaining par- 
ticle sizes for Hermsen 's data base, since they generally 
result in only 0^2 (volume -surf ace mean diameter) and are 
considered biased towards small particles. However, this 
technique has several distinct advantages. It is non- 
intrusive and theoretically can be used at any location in 
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or outside of the motor. Particle collection procedures and 
scanning electron microscope (SEM) evaluation can be used to 
determine the particle size distribution and, if this fits a 
log-probability function properly, the measured 0^2 can be 
directly related to [Ref- 8 ]. 

This investigation was part of a larger ongoing effort at 
the Naval Postgraduate School. Its scope was (a) to develop 
and calibrate a dif fractively scattered laser power appara- 
tus, (b) to apply the technique at the entrance and exit of 
the exhaust nozzle of a small rocket motor, and (c) to ob- 
tain "across nozzle" particle size data from various propel- 
lants which can be used to validate or improve upon the 
existing particle size/behavior models used in the Air Force 
Improved Solid Performance Program (SPP) . 
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II. DIFFRACTIVELY SCATTERED LIGHT METHOD 



A, INTRODUCTION 

The method has its origin in the Mie theory of light 
scattering by spherical particles. Gumprecht and Sliepcevich 
[Refs. 9 and 10] developed from this a theory describing the 
scattering properties of a polydispersion. Using this theory 
Chin, Sliepcevich, and Tribus [Ref. 11] conducted an investi- 
gation into the determination of particle size distributions 
in a polydispersion. This theory was limited to a small op- 
tical depth and there were limits on particle size and re- 
fractive index [Ref. 12]. A more general theory, applicable 
to finite optical depths, arbitrarily large particle size and 
arbitrary refractive index, was presented by Dobbins, Crocco 
and Glassman [Ref. 12]. However, their theory was applicable 
only to certain types of droplet size distributions. Roberts 
and Webb [Ref. 13] alleviated this problem and demonstrated 
the applicability of the "upper-limit distribution function" 
(ULDF) which was introduced by Mugele and Evans [Ref. 14] . 
Using the ULDF Roberts and Webb concluded that the volume- 
surface mean diameter, D^-, of a polydispersion could be 
accurately determined from the intensity of dif fractively 
scattered light without any knowledge of the general dis- 
tribution type. Dobbins and Jizmagian [Refs. 15 and 16] 
also concluded that the mean scattering cross section of 
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particles was primarily dependent on the mean diameter D ^2 / 
and only weakly dependent on the shape of the size distribu- 
tion function when the refractive index is near unity. The 
ratio of scattered light at two forward angles was shown to 
be relatively insensitive to particle refractive index and 
concentration by Hodkinson [Ref. 17]. He also demonstrated 
that this ratio yields a measure of 032 - 

With the above theories as a basis, experimental methods 
have been developed for estimating the volume -surface mean 
diameter , 0^2 ^ of a polydispersion of particles. Nejad et al 
[Ref. 18] used this method for measurement of the mean drop- 
let diameter resulting from atomization of a transverse 
liquid jet in a supersonic airstream. Mean diameter, refrac- 
tive index and volume concentration of aerosols have been 
determined by Powell et al [Ref. 19] using this technique in 
conjunction with light transmission measurements. More re- 
cently the practicality of this method for determination of 
particle sizes in the exhaust of a solid propellant rocket 
has been generally demonstrated by Karagounis [Ref. 20]. 

B. TECHNIQUE 

The theory for the scattering of light by a single di- 
electric spherical particle of size number, a, gives the fol- 
lowing expression for the radiant intensity 1(6) [Ref. 12]: 
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where D is the diameter, a is the particle size number (ttD/A) , 
\ is the wavelength of incident light, m is the refractive 
index of the scattering media, 9 is the angle measured from 
the forward direction (radians) due to an incident planar 
wave of irradiance Eq and is the Bessel function of first 
kind of order unity. 

The three terms in the bracket of Equation (1) represent 
the Fraunhofer diffraction, the optical scattering due to 
refraction of the centrally transmitted ray and the optical 
scattering due to a grazingly incident ray, in that order. 

The following restrictions apply to Equation (1) [Ref. 12]: 

(1) The incident light must be planar and monochromatic. 

(2) The forward angle 9 must be small (9 s sin 9) . 

(3) The particle size number (a = ttD/X) , and phase shift 
(p = 2a(m-l)) should be large. 

(4) The distance between particle and observer should be 
large compared to D^/X. 

(5) The particle must be non-absorbing. 

Dobbins et al [Ref. 12] found that for a polydispersion, 
when the expression for intensity of scattering is normalized 
(by dividing by the intensity of dif fractively scattered 
light in the forward direction (9 = 0)), the second and third 
terms in Equation (1) were small and could be ignored. Thus, 
the expression for the normalized integrated intensity of 
forward scattered light, 1(9), due to a polydispersion of 
large particles is given as: 



16 



I (6) 



( 2 ) 



D 



n 

^ L 



c» 



2J, (a0) 



N (D) D^^dD 



0 



/ Deo 4 

N (D)D dD 

0 



where is maximum diameter expected; and N(D) is a distri- 
bution function such that the integral of N(D) over the di- 
ameter interval is the probability of occurrence of particles 
within that interval. Equation (2) is valid when all parti- 
cles are illuminated equally or when the attenuation of the 
incident beam is slight. 

The transmission law for a polydispersion of particles 
is given by [Refs. 12, 15, and 21]: 



where Q is the average extinction coefficient; C and C are 
the volume and mass concentration of particles, respectively, 
T, turbidity; L, path length, and p, particle density. Q and 
0^2 sre defined as [Ref. 22]: 




(3) 



and by [Ref . 22 ] : 
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Dobbins [Ref. 12] reconunends that the optical depth (tL) 
be maintained below 1.5 to prevent distortion of the illumina- 
tion profile, 1(9) vs 8. 

Using the ULDF developed by Mugele and Evans [Ref. 14] , 
Roberts and Webb [Ref. 13] have concluded that a value ,of D ^2 
may be determined from the intensity of dif fractively scat- 
tered light from a polydispersion of spherical particles, for 
a wide range of size distributions. They developed a uni- 
versal illumination profile (shown in Figure 5) which is a 
mean curve showing the relation between normalized intensity 
1(6), and 9, a reduced angle of scattering equal to 
Using this curve one can determine mean particle size but not 
the distribution of sizes. 

The value of 0^2 then be determined as follows. 

(1) The optical apparatus is calibrated by using parti- 
cles of known size as a scattering media, i.e. 1(6) vs 6 
curves are generated and compared to the universal curve. 

(2) The actual particles, in the current investigation 
the exhaust of the rocket motor, are introduced and 1(9) vs 
6 is determined. 

(3) Knowing 1(9) and 9 , 9 is determined from the cali- 
bration curve and/or the universal curve. 

(4) D ^2 is determined using the relation; 
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C. APPLICABILITY 



Karagounis [Ref. 20] has demonstrated the applicability 
of this technique in estimating the size of Al^O^ particles 
in the exhaust of a rocket motor. The restrictive conditions 
for the applicability of equation (2) were satisfied in the 
present experiment in a similar manner. 

(1) A He-Ne laser with spatial filter and beam expander- 
collimator provided a planar, monochromatic light source. 

(2) The forward angle 9 was less than 6°. 

(3) Large particle size number (>1) a = itD/A was satis- 
fied (with X = .6328 ym, 0.2 ym < D < 20 ym, then 1 < a < 100) . 

(4) Large phase shift was provided (>1) . (Estimating 
m = 1.76 for AI 2 O 2 [Ref. 23] and 1 < a < 100, 1.52 < p < 152.) 

(5) The distance between the photodiode array and the 

particles was set at .8m, which is considerably greater than 

D^/X = 632 ym (for a D of 20 ym) . 

max 

(6) The final condition calling for non-absorbing parti- 
cles is not met exactly, however the absorbing coefficient 
for Al^O^ is approximately 10^ times smaller than the scat- 
tering coefficient for the wavelength, refractive index and 
size numbers considered here [Ref. 20]. Plass [Refs. 23-25] 
presents a very detailed discussion of this aspect. 

The additional requirement of optical depth, tL < 1.5 is 
satisfied by not applying this technique to situations where 
E/Eq < 0.25 (Eq. 4) . 
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III. 



EXPERIMENTAL APPARATUS 



A. LIGHT SCATTERING APPARATUS 

A schematic diagram of the final optical set-up is shown 
in Figure 1 and photographs are shown in Figure 2 . A sche- 
matic of the Data Acquisition System is shown in Figure 3. 

The instruments and optical equipment are listed in Table I. 

A complete description of the Data Acquisition System is pre- 
sented in Reference 26. The optical equipment was mounted 
on two parallel optical benches, one for the beam passing 
through the exhaust and one for the beam passing through the 
motor at the nozzle entrance. Initially a 5 mw Helium-Neon 
laser was used as the light source. This was later replaced 
by a 10 mw Argon laser, for reasons to be discussed later. 

A spatial filter/beam expander located directly in front of 
the light source produced a uniform light beam, one centi- 
meter in diameter. Initially a pinhole with a 2.2 mm diame- 
ter was located in front of the collimating apparatus. This 
pinhole was used to reduce the size of the laser beam in 
order to decrease the measurement volume. Further experi- 
mentation eventually determined that this pinhole was not 
necessary. The intensity of the beam was decreased using a 
1 percent neutral density filter to prevent saturation of 
the linear photodiode arrays. The beam was split using a 
50/50 cube beam splitter. One beam continued through the 
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exhaust, the other traveled to a 90 degree prism where it was 
turned parallel to the first beam and directed through the 
motor at the nozzle entrance. Passage of the beam through 
the motor was afforded by mounting circular glass windows in 
the walls of the motor on either side of nozzle entrance 
area. Both beams were passed through identical condensing 
lenses (5 cm in diameter with focal lengths of 500 mm) and 
focused on linear photodiode arrays. The photodiode array 
assemblies each consisted of a linear photodiode array, a 
"mother board" circuit card and an array board. Each as- 
sembly was mounted in a lighttight black box fitted with a 
narrow band pass filter. Each box was mounted on two trans- 
lation stages to provide vertical and horizontal degrees of 
freedom. The linear photodiode array was a self-scanning 
photodiode array consisting of 1024 silicon photodiodes 
mounted in a vertical row on 25 ym centers. The aperture 
width was 26 ym. The circuits provided an integrated, 
sampled-and-held output. The "mother board" contained most 
of the circuitry, including the driver/amplif ier , clock, 
start, and end of scan functions. The array board contained 
only the circuitry which must be located close to the array 
itself. The output signal was passed through a variable low 
pass filter (for reasons discussed later) before transmission 
to the data acquisition system. The reading of the array was 
accomplished in 36 milliseconds, which corresponds to a rate 
of 30 kHz. 
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B . ROCKET MOTOR 



The motor was cylindrical, stainless steel with a copper 
nozzle. The chamber was two inches in diameter and two 
inches deep (Figure 4a) . The nozzle throat diameter was .245 
inches (for a chamber pressure of 760 psi) . An internally 
burning, six-pointed star (Figure 4b) was selected for the 
grain design in an effort to achieve a period of steady state 
pressure in which to take data. The grains were one inch 
long and two inches in diameter with a web thickness of 0.5 
inches. Each grain was cut out manually from a one inch 
thick slab of propellant using frabicated shaped cutters. 

The grains were prepared for ignition by applying a thin 
black powder and glue mixture to the inner surface of the 
grain. Electrically ignited BKNO^ was discharged from the 
head end of the motor onto the inner surface of the grain, 
thus igniting the black powder, which in turn ignited the 
propellant. A cylindrical stainless steel tube (8" O.D.) was 
mounted at the exhaust of the motor to collect the samples of 
exhaust particles for SEM evaluation. 

Two circular windows were mounted in the walls of the 
motor on either side of the nozzle entrance area. The win- 
dows were recessed from the motor chamber. To keep the win- 
dows clean, nitrogen was discharged into this recessed area 
and evenly diffused through a sintered metal filter (Figure 
4c) . 
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IV. EXPERIMENTAL PROCEDURE 



The apparatus was established in an open test cell of the 
Combustion Laboratory. The output signal of the system con- 
tained a quantity of high frequency noise. After numerous 
unsuccessful attempts to locate and eliminate the source of 
this noise, it was decided to use a variable low pass filter 
to eliminate the noise from the output signal. Subsequent 
examination of the filter effects during the calibration runs 
determined that the filter did not adversely affect the ac- 
curacy of the data readings. After a collimated and filtered 
beam was acquired, the focus was set 1.52 mm above the first 
diode of the array. Care was taken at all times to avoid 
focusing the beam directly on the diodes to prevent damage to 
the array. It was also found necessary to slightly cant the 
focusing lens to prevent interference from reflected light. 

Next, the accuracy of the Dif fractively Scattered Light 
Method (DSLM) for obtaining the mean diameter of particles 
was investigated. A comparison was made between experimen- 
tally determined normalized intensity (1(9)) versus non- 
dimensional angle (^) profiles and the "universal" theoretical 
profile (Figure 5) . The various profiles were first obtained 
using spherical glass beads of different size ranges. The 
individual profiles for the three size ranges of beads (1-37 
ym, 37-44 pm and 53-63 pm) are shown in Figures 6-8. 



23 



The profiles for each size range were determined by sus- 
pending the particles in water contained in a Plexiglas box 
(2 . 5” X 3 . 0 ” X 2 . 5" ) . Suspension was maintained using a mag- 
netic stirrer. The first step in this calibration procedure 
was to obtain the zero-scattering profile, that is, the light 
intensity on each diode after passage through the stirred 
water, before particles were introduced. Next the particles 
were introduced and the intensity reading at each diode was 
taken. Examples of intensity versus diode at 1 and 3 kHz 
filter settings are shown in Figures 9 and 10, respectively. 
The analog to digital conversion was accomplished by the data 
acquisition system. (A complete discussion and listing of 
the data acquisition and the reduction procedures discussed 
below are presented separately by Hansen [Ref. 26].) 

After obtaining the two scattering profiles the normal- 
ized intensity 1(0) versus angle theta (6) profile was de- 
termined according to the relation: 



1 ( 0 ) 



In - In V- - 

9/P 9/0 _ 0,p 0,0 

I , - I , V ^ - V T 

cl,p cl,o cl,p cl,o 



( 8 ) 



where or are the projected light intensity and re- 

lated voltage at the centerline, and Ig and Vg are the light 
intensity and voltage at an angle, 0 (Figure 11) . "o” and "p" 

symbolize no particles and with particles, respectively. The 
first step in determining the normalized intensity was to 
calculate the nximerator, /^Vg , of equation (8) by subtracting 
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the two voltages at corresponding diode numbers. The maximum 
difference (lower limit on diode number) computed was taken 
as the point where the dif fractively scattered light inten- 
sity was dominant over the transmitted light intensity. The 
computer determined the diode number where this occurred. An 
upper limit on diode number was taken near the point where 

Va - Va first became zero (Figure 10) . A least squares 
y y p y ^ o 

linear fit was then made to the AVg vs diode number data and 
extrapolated to obtain ^ ^ (the forward scattered 

light at 9 = 0) . 1(6) vs 9 could then be calculated (equa- 

tion (8) ) . 

An example of the relationship of 1(9) versus 9 thus pro- 
duced is shown in Figure 12. At this point the relationship 
of 1(9) versus 9 is determined using the relationship: 



A 

where 0^2 is the mean diameter of the particle range being 
used in the calibration. 1(9) vs 9 can then be plotted (Fig- 
ure 13) . An nth order polynomial curve fit was then applied 
to the data over a user specified range. The order of the 
polynomial (also user specified) may be varied from 1 to 10 
depending on the appearance of the data. The range corres- 
ponded to the range of the diodes selected for determination 
of the centerline values as mentioned above. An example of 
the result of this process is shown in Figure 13. The 
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resulting curve is then plotted on a semi-log scale for com- 
parison with the universal curve (Figures 5-8) . 

The results shown in Figures 5-8 for glass spheres gener- 
ally agree with those of Karagounis [Ref. 20] in that larger 
deviations from the universal curve occur as 9 increases (be- 
yond ¥ = 2.5) . However, unlike the data of Karagounis, 
agreement with the universal curve was not dependent on 
particle mean diameter. 

Several possible causes for this lack of agreement were 
investigated. First the effect of filter setting was checked. 
Keeping all other variables constant, data were taken with 
1 kHz and 3 kHz filter settings. The resultant profiles 
shown in Figure 14 are almost identical. Next the effect of 
horizontal displacement of the array was studied. Again 
keeping all other variables constant, two runs were made. 

The first had the focal point centered on the axis of the 
array. In the second, the focal point was displaced to one 
side just far enough to minimize non-uniform intensity varia- 
tions. Again the results showed little difference from one 
another (Figure 15) . A narrower particle suspension box was 
also used to decrease the depth of the field, but it likewise 
showed little change from the profile using the wider box. 
Signal-to-noise ratio and dark current were two factors which 
could not be completely eliminated and may contribute to the 
lack of agreement at higher 9. 
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Using the same procedure as discussed above, 1(0) versus 
'd profiles were prepared for four size ranges of non-spherical 
AI 2 O 2 powder. The results are shown in Figure 16. The re- 
sults indicate reasonable agreement with the universal curve. 
The experimental profiles had greater slopes than the theo- 
retical profile as they did for the spherical glass beads. 
However, for the AI 2 O 2 the mean size affected the profile. 

The particle size measuring technique was then applied 
to actual motor firings. Before each motor firing calibra- 
tion of the recording equipment was conducted. Proper posi- 
tioning of the laser beam on the array was verified. Then 
the zero scattering profile was recorded. The motor was 
fired and the scattering profile was recorded automatically 
at a specified pressure. Chamber pressure was recorded on a 
visicorder oscillograph. In the same manner previously de- 
scribed the normalized intensity (1(9)) versus angle (9) was 
determined. 0^2 values were then calculated over the re- 
corded 9 range, using both calibration curves and the uni- 
versal curve, from the relation: 



D 
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9tt 



( 10 ) 



Residue from the motor exhaust was captured in a stain- 
less steel collection tube. Some of it was dissolved in 
acetone, subjected to an ultrasonic vibrator, allowed to 
settle and then dried in a vacuum oven. This "cleaned" 



27 



sample and an "as collected" sample were examined with an 
SEM. Photographs are shown in Figures 17 and 18. 
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V. RESULTS AND DISCUSSION 



The original goals of this investigation were (a) to 
develop an apparatus to determine D ^2 from measurements of 
the dif fractively scattered laser power spectra using rapid 
scans of linear photodiode arrays, (b) to apply the apparatus 
at the nozzle exit and through the motor at the exhaust noz- 
zle entrance of a solid propellant rocket motor, and (c) to 
obtain "across-nozzle" particle data for comparison with the 
existing models for two-phase flow phenomena in exhaust noz- 
zles. The first two goals were met, however, due to experi- 
mental difficulties the last goal was only partially achieved. 

The apparatus developed for this investigation has been 
discussed above. Considerable difficulty was encountered 
during the initial set-up of the computer controlled data 
acquisition and data reduction procedures. The resulting 
experimental procedures have also been discussed above (Sec- 
tion IV) . When the apparatus was applied to actual motor 
firings several additional difficulties occurred. These 
will be discussed below. 

Keeping the windows clear without affecting either motor 
flow or the combustion processes presented the greatest 
challenge. The final solution (Figure 4c) consisted of a 
sintered metal filter to evenly diffuse the nitrogen into 
the window purge area and a converging nozzle to slightly 
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accelerate the nitrogen flow into the motor. This system 
worked very well at all pressures investigated. 

The light generated during the combustion process within 
the motor (but not at the nozzle exit) was found to produce 
significant radiation at the wavelength of the He-Ne laser 
(.6328 ym) initially used as the light source. This pro- 
hibited the measurement of any scattered light from the par- 
ticles. The light source was changed to an argon laser 
(.4880 ym) and the problem was alleviated. 

The location of the laser beam in the exhaust jet was 
found to be critical. If directed too far to the rear of 
the nozzle exit, it either missed the jet entirely due to 
the deflection of the jet or the particle number density was 
too low. Both problems resulted in a particle scattering 
profile that was nearly identical to the zero particle scat- 
tering profile. With the current motor the beam was located 
three inches to the rear of the nozzle exit. 

The tests conducted are summarized in Table II . The 
first series of tests was conducted at 775 psia. There 
were two significant results: (1) there was insufficient 

light through the motor to get a scattering profile; (2) the 
measured 0^2 the exhaust was large (30-50 ym) . The photo- 
graphs in Figures 17 and 18 show the apparent reason for this 
large 0^2 to be the presence of a small number of large ir- 
regularly shaped agglomerates, which were not A1 or AI 2 O 2 
(typically J?-5 ym spheres) . These agglomerates were most 
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probably inhibitor or binder char. Their influence was sig- 
nificant, as even a very few large particles greatly increase 
the measured 032 * sensitivity of the technique to ir- 

regularly shaped particles has already been demonstrated 
(Figure 16). Repeated tests resulted in different 0^2 
the nozzle exhaust, apparently from different sizes of char 
agglomerates at the time the diode scan was made. 

The tests made at approximately 500 psia (runs 3, 4 and 
5, Table II) were made to determine, (1) if the size of the 
exhaust char agglomerates were sensitive to combustion pres- 
sure, and ( 2 ) if the lower chamber pressure would permit 
measurement of scattered light through the motor. In addi- 
tion, repeated runs were made to acquire the diode data at 
different times during the run to determine whether the data 
was significantly time dependent. It was found that data 
could be obtained through the motor. Comparisons with the 
universal curve and the small glass bead calibration curve 
(Figure 6 ) yielded reasonable values for 0^2 of 16 and 12 ym, 
respectively. It was further determined that 0^2 measured at 
the nozzle exit was sensitive to burn time, as shown by the 
increasing values of 0^3 with increased time. And finally, 
as shown in Figure 17, the spherical Al/Al^O^ particles were 
still in the ^5 to 5 ym size range. Therefore, the change in 
measured 0^2 was apparently due to changes in char agglomer- 
ate size during the run. 
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In run 6 an end burning grain was used in an attempt to 
reduce or eliminate the large char agglomerates. The results 
were: (1) a significant reduction in measured 0^2 (however, 

it was still larger than the AI/AI2O2 spheres) , and ( 2 ) the 
larger optical path with particles (increased optical den- 
sity) precluded the measurement of any scattering data. 
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VI. CONCLUSIONS 



The technique developed for determining the change in 
0^2 across a solid propellant rocket motor exhaust nozzle 
appears to work well. Additional effort is required to find 
a grain configuration and/or an inhibitor that minimizes the 
size of, or eliminates entirely, the char agglomerates. 

This is required in order to be able to obtain "across noz- 
zle" data during a single test. This would allow realization 
of the final goal, that is, to obtain data for the validation 
of the two-phase flow models. 
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TABLE I 



EXPERIMENTAL APPARATUS 

1. Spectra Physics Model 147, Helium-Neon Laser. Power 
5 mW, Random Polarization. 

2. Power Supply for above, Model 247, 30 watts. 

3. Laser Beam Expander with spatial filtering, 22 mm aper- 
ture, Oriel, Model 1526. 

4. Plano-convex lens, 5 cm diameter. 50 cm focal length. (2) 

5. Neutral density filter, FNG-043, 1%, 2.0 O.D. 

6. Cube Beam-Splitter, 50/50. 

7. 90° Prism. 

8. Linear Photodiode Array, EG&G Reticon, G series, 1024 
elements. (2) 

9. Mother Board Circuit Card, EG&G Reticon, RCIOOB. (2) 

10. Array Circuit Card, EG&G Reticon, RC106. (2) 

11. Narrow Band Pass Filter. (2) 

12. DC Power Supply for Photodiode Assembly. 

13. Variable Low Pass Filter, Krohn-Hite, model 3343. 
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Propellant contained 2% by weight of A1 (40 pm) 
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Figure 1: Schematic Diagram of Dif f ractively Scattered Light Apparatus 








Figure 2: Photographs of Light Scattering Apparatus 



37 










o 

q: 

K 



O 

o 



UJ 

ac 

I 



o 

QC 

< 

X 



38 



Figure 3: Schematic of Data Acquisition System 
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Figure 4: Motor Components 
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Figure 5: 1(0) vs 0, Normalized Intensity Profiles, Spherical Glass Beads 
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Figure 6: 1(6) vs 9; Spherical Glass Beads, 
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Figure 7: 1(0) vs 0; Spherical Glass Beads, 37-44 
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Figure 8: 1(6) vs 0; Spherical Glass Beads, 53-63 
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Figure 9: Voltage vs Diode--l kHz Filter 
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Figure 10: Voltage vs Diode — 3 kHz Filter 
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Figure 11: Light Scattering Geometry 
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Figure 12: 1(0) vs 0--Example 
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Figure 13: 1(6) vs 0--Example 
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Figure 14: 1(0) vs 0 — 1 and 3 )cHz Frequency Filter Settings 
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Figure 15: 1(0) vs 6 — Horizontal Displacement of Array Comparison 
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Figure 17: 



SEM Photograph — Clean Sample 



(Typical ) 
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Figure 18: SEM Photograph — As Collected (Typical) 
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